Emission angle dependent pion interferometry at RHIC and beyond 



(N 
O 

o 

(N 



(N ' 

> : 

00 . 

r» ■ 

(N ■ 

: 

o . 

(N ■ 

o ; 

^ : 

a: 
i 

Q . 
^ ■ 



X 



Ulrich Heinz 1 ' 2 and Peter F. Kolb 2,3 

'Physics Department, The Ohio State University, Columbus, OH 43210 
Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030 
Department of Physics and Astronomy, State University of New York, Stony Brook, NY 1 1794-3800 

(Dated: July 10, 2002 (revised version)) 

We use hydrodynamics to generate freeze-out configurations for non-central heavy-ion collisions at 
present and future collider energies. Such collisions are known to produce strong elliptic flow. The 
accompanying space-time structure of the source at freeze-out is analyzed using pion interferometry. 
Between RHIC and LHC energies the source deformation in the transverse plane changes sign. This 
leaves characteristic signatures in the emission angle dependence of the HBT radii. 
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Heavy-ion collision data from the Relativistic Heavy 
Ion Collider (RHIC) jf) have produced strong indications 
for a large degree of thermalization in these collisions. 
Only dynamical models with very intense reinteractions 
among the produced particles |2L pi kL in particular the 
hydrodynamic model ]|, |[ [j], ]s|, UllOl, pLl}| , have been 
able to reproduce the large observed elliptic flow signals 
U) , [l3[ [l4| [ll| together with the measured systemat- 
ics of the shapes of the single-particle hadron spectra 
]Tg| . |l7| ] which require strong radial flow jL8| . When com- 
pared to the data, all models seem to fail which are not 
able to achieve almost complete local thermalization very 
quickly, at most 1 fm/c after nuclear impact |n| S3). At 
this early time the energy density in the reaction zone 
is still far above 1 GeV/fm 3 [Q, the critical value for 
color deconfinemcnt, so that the thermalized state formed 
early in the collision has most likely been a quark-gluon 
plasma [ pp[ . 

Even though the microscopic mechanisms causing the 
fast thermalization are presently not understood, these 
observations provide strong phenomenological support 
for the hydrodynamical model as a dynamical description 
for the space-time evolution of the reaction zone. We will 
therefore use it here to study additional aspects of the 
collision. We will concentrate in particular on the ques- 
tion how the strong elliptic flow affects the spatial matter 
distribution at "freeze-out" when the observed hadrons 
decouple from the fireball. Elliptic flow requires an initial 
spatial deformation of the nuclear overlap region in the 
transverse plane (i.e. either a non-zero impact parameter 
or collisions between deformed nuclei) since it is gener- 
ated by anisotropies in the pressure gradients || |l| . As 
it develops, the matter begins to expand more rapidly in 
the "short" direction (which for non-central collisions be- 
tween spherical nuclei lies in the reaction plane). Hydro- 
dynamics predicts that, if the initial energy density is 
high enough and thus the time until freeze-out sufficiently 
long, this initial out-of-plane deformation changes sign, 
leading to an "in-plane extended source" (IPES) at de- 
coupling. We show that at RHIC energies this is not yet 
the case p2| , but that it will eventually happen at higher 



collision energies (LHC or beyond) . In this paper we dis- 
cuss how this transition can be probed by two-particle 
interferometry, in particular by studying the emission an- 
gle dependence relative to the reaction plane of the size 
parameters ("HBT radii") extracted from Bose-Einstein 
correlation functions [||[ |24|, [2(| ^7j. We show that 
the spatial and momentum anisotropies of the phase- 
space distribution at freeze-out (the emission function 
S(x, K)) generate characteristic azimuthal oscillations of 
the HBT radii whose phases and amplitudes reflect this 
transition. The way in which they reflect it turns out to 
be non-trivial and will be discussed in some detail. 

The hydrodynamic model and equation of state used 
in our calculations are described elsewhere @|. We here 
study semicentral Au+Au collisions at impact param- 
eter 6 = 7fm. b is taken to point in x direction, z is 
the beam axis such that (x, z) span the reaction plane, 
and y points perpendicular to it. We assume boost- 
invariant longitudinal expansion and concentrate on the 
dynamics in the transverse (x,y)-plane at z = where 
the nuclei collide. We use two different sets of ini- 
tial conditions. The first set (labelled "RHIC1") corre- 
sponds to optimized values for central Au+Au collisions 
at y/s = 130 A GeV, which also give a good description 
of the measured hadron spectra and elliptic flow in non- 
central collisions @, [H| ^ [D| |l6| [f| out to transverse 
momenta of about 2GeV/c and impact parameters of 
about 10 fm ||, |20| . The initial entropy density distri- 
bution is parametrized by a superposition of two compo- 
nents which scale with the number of wounded nucleons 
(with 75% weight) and binary nucleon-nucleon collisions 
(with 25% weight), respectively (^8|. The corresponding 
peak value of the initial temperature, at the starting time 
to = 0.6 fm/c of the hydrodynamic expansion, in central 
(6 = 0) collisions is To = 340 Me V, scaled down in semi- 
central collisions with the Glauber prescription described 
in pj|. Freeze-out is implemented with the Cooper- Frye 
algorithm p3| along a surface of constant temperature 
T f = 130MeV. 

The second set of initial conditions (labelled "IPES" ) 
uses a much higher initial temperature Tq — 2 GeV (at 
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b = 0), and the hydrodynamic evolution is started at 
To = 0.1fm/c and stopped when a freeze-out tempera- 
ture Tf = 100MeV has been reached. Such a high ini- 
tial temperature can probably not even be achieved at 
the LHC, and the strong transverse flow generated in 
this case probably causes the system to decouple already 
much closer to the hadronization temperature of (in our 
case) Th a d = 164 MeV. We force the system to start at 
such a high temperature and freeze out so low in order 
to give it enough time to convert the initial out-of-plane 
deformation into a final in-plane deformation (which is 
the phenomenon we want to study) . Unfortunately, over 
large times and due to the strong radial flow numeri- 
cal instabilities in the hydrodynamic code accumulate, 
and to keep these at a minimum we use (for IPES ini- 
tial conditions only) an equation of state (EOS) without 
a sudden phase transition. Arguing that with these ini- 
tial conditions the system anyway spends almost all of 
its time in the quark-gluon plasma phase, we simply use 
P=^e, corresponding to an ideal gas of massless quarks 
and gluons. After stopping the evolution at Tf, we as- 
sume that pions are emitted with this freeze-out tem- 
perature and the hydrodynamic flow established on the 
corresponding freeze-out surface. This implies a change 
of EOS at freeze-out and is thus not entirely dynamically 
consistent; however, since modifications of the dynamics 
resulting from a proper transition to a hadron resonance 
gas EOS at hadronization are expected to have at most 
minor effects on the source anisotropy at such a late time, 
this dynamical inconsistency should not qualitatively af- 
fect our results. 

The pion emission function S(x, K) for pions of a given 
charge is computed from the hydrodynamic output by the 
Cooper-Frye prescription p9[ 



S{x,K) 



(2tt) 3 ,/ s ex.p[K-u{x')/T t ] - 1 



(1) 



This includes only directly emitted pions. While reso- 
nance decay pions strongly affect the shape of the single- 
particle spectrum, the same is not true for the HBT radii 
extracted from the two-pion correlator Q. They do, 
however, generate a non-Gaussian spatial tail in the emis- 
sion function ]9|, |30| which significantly increases its spa- 
tial widths ("homogeneity lengths", see p^j), thereby de- 
stroying the direct correspondence between the latter 
and the measured HBT radii which exists for Gaussian 
sources ^6|. Since we want to preserve this correspon- 
dence as much as possible, we exclude decay pions when 
calculating the homogeneity lengths. We have checked 
that the remaining differences between the HBT radii 
calculated from the homogeneity lengths (see below) and 
those extracted from a Gaussian fit to the correlation 
function (calculated in the standard way Q by Fourier 
transforming the emission function) are small enough to 
not affect the qualitative features of our results. 

In equation ([!]), £ denotes the freeze-out surface U(x) 
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FIG. 1: (a) Cuts through the freeze-out surface U(r) at z~0 
along and perpendicular to the reaction plane, (b) Contour 
plots in the transverse plane of the time-, z-, and momentum- 
integrated emission function d 2 N/d 2 r for RHIC1 and IPES 
initial conditions (see text), (c) Cuts through diagram (b) 
along and perpendicular to the reaction plane, for RHIC1 
and IPES initial conditions. 



of constant temperature Tf , and u^(x') is the flow veloc- 
ity on £. Figure [l^, shows cuts through £ along (y — 0) 
and perpendicular to the reaction plane (x = 0) . Initially 
the source is extended out-of-plane (larger in y than in 
x direction), but it then expands more rapidly into the 
x-direction, becoming in-plane elongated at later times. 
For RHIC1 initial conditions this only happens after most 
of the matter has already decoupled; as a consequence the 
time-integrated source, shown in Fig. [j]b, is still longer in 
y than in x direction in the RHIC1 case. For IPES initial 
conditions the deformation changes sign before most par- 
ticles decouple, and the time-integrated source appears 
in-plane-extended (see again Fig. [l]b). Also, it is much 
larger due to the much higher initial energy density and 
longer lifetime. 

Figure |l|c shows cuts along and perpendicular to 
the reaction plane through the density contour plots 
d 2 N/d 2 r= J(d 3 K/E K )dzdtS(x,K) of Fig. gb. Pion 
emission is seen to be strongly surface peaked, in par- 
ticular at RHIC1 where the freeze-out radius is almost 
constant for a long time. This "opacity" is weaker both 
at lower collision energies (where the freeze-out surface 
shrinks to zero continuously |]33|) and at higher energies, 
due to larger temporal variations of the freeze-out radius. 

Figures ^ and || show the spatial distributions of pi- 
ons emitted with fixed momentum. Shown are density 
contours of J dz dt S(x, K) in the transverse plane (x, y) 
for pions with rapidity Y — and fixed K±, for three 
emission angles $ = 0, 45° and 90° relative to the reac- 
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FIG. 2: Contours of constant emission density in the trans- 
verse plane for RHIC1 initial conditions (see text). The thick 
line indicates the largest transverse extension of the freeze- 
out hypersurface (see Figs. |l|a,b). The four panels show 
emission regions for midrapidity pions (Y = 0) with K± = 
and, for three emission angles indicated by arrows, with 
K ± =0.5GeV. 
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FIG. 3: Same as Fig. ^, but for IPES initial conditions. 



tion plane. (In Fig. || we replaced $ = 45° for clarity 
by the equivalent angle <& = 225°.) Particles with van- 
ishing transverse momentum K± are seen to be emitted 
from almost the entire interior of the "bathtub" shown 
in Figs. |b,c; for RHIC1 (IPES) initial conditions this 
region is elongated out-of-plane (in-plane). For slow pi- 
ons the source thus looks transparent. Pions with suf- 
ficiently large transverse momenta are emitted from rel- 
atively thin regions close to the rim of the "bathtub" 
where the flow velocity is largest and points into the di- 
rection of the emitted pions. For fast pions the source 
thus looks opaque |34[ [35l [3(|. Their emission regions 
rotate with the emission angle, constrained by the shape 
of the "bathtub" even though they are much smaller. 
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FIG. 4: Azimuthal oscillations of the HBT radii at Y — for 
6=7fm Au+Au collisions at yfs = 130 A GeV (RHIC1), for 
four values of the transverse momentum K i as indicated. 



Again we see that the opaqueness is stronger at RHIC 
than at higher energies. It is also seen to be anticorre- 
lated with the curvature of the "bathtub" wall: where 
the curvature is large, the source is squeezed less tightly 
to the wall than where it is smaller. Fig. ||b shows that 
for IPES this geometric effect wins over the anisotropic 
flow effects which push the emission region more strongly 
towards the wall in x than in y direction; for RHIC1 both 
effects act together. 

We now proceed to compute the HBT radii from the 
widths of the if-dependent emission regions [^6). The 
indices o (for "outward") and s (for "sideward") indi- 
cate the directions parallel and perpendicular to the emis- 
sion vector K±; I indicates the longitudinal z direction 
p6| . Given the spatial correlation tensor S^ u — {x^x v ) — 
{x IM ){x v ), which describes the (if- dependent) widths of 
the emission region in space and time, the HBT radii can 
be calculated at midrapidity (Y = 0) from the relations 
MM 



R s 

K 

1X1 OS 



S x 



$xx $yy 



^cos(2$)-S , a; j / sm(2$), 
cos(2$) + S xy sm(2&) 



-2p x (Stx cos $+S ty sin $) + (3 2 ± St. 



S xy cos(2$) 



sin(2$) 



-P±(S tx sin ®-S ty cos $) 



(2) 



Note that R% and i?^ s receive purely geometric (first 
lines) and mixed space-time correlation contributions 
(second lines of the equations). Since we want to use 
pion interferometry to obtain information on the geomet- 
ric deformation of the source at freeze-out, these must be 
discussed separately (see Fig. || below). 

In addition to the explicit $ dependence exhibited in 
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FIG. 6: The tilt angle of the major axes of the emission 
region in the transverse plane relative to the (o, s) system, for 
different values of K i . 



FIG. 5: Same as Fig. ^, but for IPES initial conditions. For 
i?o and _Rq S the geometric contributions are shown separately 
as thin circled lines (see text). 



Eqs. (g), which arises from the angle <& between the (o, s) 
and (x, y) coordinate systems p4| , the spatial correlation 
tensor S^„ depends implicitly on $ through the emission 
function 5(x; Y, $). The combined $ dependence 
of the radius parameters is shown in Figs. || and for 
RHIC1 and IPES initial conditions, respectively. We first 
note that Rf does not show any interesting azimuthal 
oscillations and requires no further discussion. The same 
turns out to be true for the emission duration Su (the 
last term contributing to R 2 in Eqs. (j|)). Interesting 
temporal contributions to the oscillation patterns of R 2 
and R 2 S thus only involve the terms S tx and S ty which 
correlate the freeze-out positions with time. 

Let us first discuss the geometric contributions. From 
Figs. ^ and || we already know that pions with K± = 
are emitted from almost the entire fireball and thus probe 
the different sign of the spatial deformation of the total 
(momentum- integrated) RHIC1 and IPES sources shown 
in Fig. [I]. This is reflected by the opposite sign of the os- 
cillation amplitudes of R 2 and of the geometric contribu- 
tion to R 2 in Figs. [| and [| at K± = 0. (For clarity we do 
not show separately the purely geometric contributions 
to R 2 and R 2 S in Fig. ^ because at RHIC1 they oscil- 
late in the same way as the total radii, in contrast to the 
IPES case, Fig [|) At higher if^-values the oscillations 
in Fig. H for R 2 S change sign, but those of the geometric 
contribution to R 2 do not. This reflects an intricate in- 
terplay between geometric and flow effects, including the 
already mentioned weaker opacity for IPES. 

The geometric contribution to the azimuthal oscilla- 
tion of R 2 S has an interesting intuitive interpretation: 
Figs. [2] and || show that for K± 7^ the emission region 
rotates together with the emission direction $, but not 
quite in phase. For each direction $ we can diagonalizc 
the spatial correlation tensor in the transverse plane. Let 
us denote the major axes by X and Y, with X pointing 



approximately outward and Y approximately parallel to 
the fireball surface. It is then easy to show that 



^(geom) = U Sxx _ 5yy)sin (2e), 



(3) 



where 9 is the tilt angle between the outward-pointing 
major axis X and the emission direction K± (resp. the 
direction of x Q ). Positive values of 9 correspond to coun- 
terclockwise rotation. Figs. || and ^ show that for K± ^ 
the deformation h(Sxx — Syy) of the emission region in 
its major axis frame is negative for RHIC1 but positive 
for IPES. For RHIC1 R 2 S is always positive in the first 
quadrant, while for IPES it is mostly negative. Accord- 
ing to Eq. (||) the tilt angle 9 is then negative in both 
cases, i.e. in the first quadrant the major axes (X, Y) of 
the emission region are always rotated clockwise against 
the (o, s) axes. This is confirmed by visual inspection 
of Figs. H and [| Figure ^| shows the tilt angle 9 as a 
function of the emission direction <E>, for various values of 
K± . The case K± = is special since the corresponding 
source is ^-independent. Then the tilt angle is exactly 
6 = —$ (mod 7r). 

Last but not least we consider the contributions from 
the space-time correlations Stx and St y to the oscillation 
amplitudes of R 2 Q and R 2 S . They are multiplied by 
(see Eqs. (0)) and thus grow with increasing K±. They 
reflect the change (caused by the elliptic flow) of the 
spatial source deformation with increasing time. Conse- 
quently their sign does not change from RHIC1 to IPES, 
contrary to the spatial contributions to the oscillations. 
For RHIC1 they reinforce the geometric oscillations of R 2 
and R 2 S , but Fig. || shows that in the IPES case they act 
against them and are responsible for changing the sign 
of the oscillations for K± £ 100 — 150 MeV. We interpret 
this change of sign of the oscillations with increasing K± 
as a possible experimental signature for the spatial man- 
ifestation of elliptic flow, i.e. the fact that the source 
expands faster in the x than in the y direction. 

Summary. We have presented a detailed analysis of the 
space-time geometry of non-central relativistic heavy ion 
collisions within an ideal hydrodynamic framework. Al- 
though the idealized hydrodynamic picture with Cooper- 
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Frye freeze-out does not Q give a quantitative descrip- 
tion of the experimental HBT radii for central Au+Au 
collisions at RHIC |37| and should eventually be re- 
placed by a more realistic approach to kinetic freeze-out 
we expect the qualitative geometric features of 
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non-central collisions discussed here and their manifesta- 
tion in the emission angle dependence of the HBT radii 
to be insensitive to details of how the particles decouple 
from the fireball. We have shown that at high enough col- 
lision energies the anisotropic flow will lead to a source 
which in a transverse cut is elongated along the reaction 
plane, orthogonal to the initial out-of-plane elongation 
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of nliclear overlap region. As an unambiguous signature 



for this change ot spatial dclormation we identified a sign 
change of the oscillation amplitude as a function of emis- 
sion angle of the sideward HBT radius R s at small trans- 
verse momentum K±. Along with this we point out ad- 
ditional changes in the oscillation patterns of R and R^ s 
which are caused by a combination of geometric spatial 
and temporal aspects of the source at freeze-out. A sign 
change as a function of increasing K± in the oscillations 
of at high collision energies has been interpreted as 
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measurable evidence for the faster growth of the source 



along the reaction plane than perpendicular to it, due to 
elliptic flow. According to our calculations the effective 
source at RHIC energies should still be elongated out of 
the reaction plane at freeze-out. The corresponding pre- 
dicted oscillation pattern of the HBT observables shown 
in Fig. ^|is in qualitative agreement with preliminary data 
from RHIC p8| ]. Qualitatively similar conclusions as to 
the sign of the source deformation at RHIC have been 
drawn by Lisa and Wells p9| , using a hydrodynamically 
motivated exploding source model to fit the spectra fll] 



and HBT radii |37 
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